Quantum-defect theory of resonant charge exchange 



(N 

o 

(N 



Ming Li 12 'E and Bo Gao 1 ' 2 ^ 

1 Department of Physics and Astronomy, University of Toledo, Mailstop 111, Toledo, Ohio 43606, USA 
2 State Key Laboratory of Low- Dimensional Quantum Physics, 
Department of Physics, Tsinghua University, Beijing 100084, China 
(Dated: May 22, 2012) 

We apply the quantum-defect theory for —1/R 4 potential to study the resonant charge exchange 
process. We show that by taking advantage of the partial-wave-insensitive nature of the formulation, 
resonant charge exchange of the type of 1 S+ 2 S can be accurately described over a wide range of 
energies using only three parameters, such as the gerade and the ungerade s wave scattering lengths, 
and the atomic polarizability, even at energies where many partial waves contribute to the cross 
sections. The parameters can be determined experimentally, without having to rely on accurate 
potential energy surfaces, of which few exist for ion-atom systems. The theory further relates 
ultracold interactions to interactions at much higher temperatures. 
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PACS numbers: 34.10.+x,34.70.+e,34.50.-s,34.50.Cx 



I. INTRODUCTION 

With the recent emergence of cold-ion experiments, ei- 
ther with trapped ions [IHl] or in the context of cold 
plasmas @, there is a growing interest in ion-atom in- 
teractions at cold temperatures [7-9], including, in par- 
ticular, the resonant charge exchange process (see, e.g., 
Refs. [|M1), such as 



Na+ + Na 



Na + Na 4 



Despite being one of the simplest reactive processes that 
has been a subject of study for a long time 0, quan- 
titative understanding of resonant charge exchange re- 
mains difficult, especially at cold temperatures. This dif- 
ficulty stems from its sensitive dependence on the poten- 
tial energy surfaces (PES), a common difficulty shared by 
all heavy particle interactions at cold temperatures (see, 
e -g-j uM)i including not onl y io n-atom interaction, but 
also atom-atom interactions [15] . and chemical reactions 
whenever the Langevin assumption breaks down [l6l . [l7| . 

In the case of atom-atom interaction, this difficulty 
has only been overcome by incorporating a substantial 
amount of spectroscopic data, especially data close to 
the dissociation limit, to fine tune the PES (see, e.g., 
Ref. [HI, HH). Without such fine tuning, no ah initio 
PES for alkali-metal systems has been sufficiently accu- 
rate to predict the scattering length and other scatter- 
ing characteristics around the threshold. The availability 
of such data, however, is limited mostly to alkali-metal 
atoms and a few other species that can be cooled. For 
ion-atom systems, with a few exceptions that came close 
[20j , no such data are yet available, though recent efforts 
on the trapping and cooling of molecular ions (see, e.g., 
Refs. |2lT[24| ) show considerable future promise. This 
status on the ion- atom PES is such that at the moment, 
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with the possible exception of H + +H and its isotopic 
variations [IJ US Hfl; n0 other predictions for cold or 
ultracold ion-atom processes, including resonant charge 
exchange, can yet be trusted before experimental verifi- 
cation. 

We present here a quantum-defect theory (QDT), not 
only as a general approach to ion-atom interactions, but 
also as one specific method of dealing with this difficulty 
of sensitive dependence on PES. It is an initial applica- 
tion of the QDT for -1/i? 4 potential 0,113, SH, as formu- 
lated in Ref. 8] , to the resonant charge exchange process. 
We show that by taking advantage of the partial-wave- 
insensitive nature of the QDT formulation 0, [2j| [3fj |. 
resonant charge exchange of the type of 1 S+ 2 S, applica- 
ble to Group IA (alkali), Group IIA (alkaline earth), and 
helium atoms in their ground states, can be accurately 
described over a wide range of energies using only three 
parameters, such as the gerade and the ungerade s wave 
scattering lengths, and the atomic polarizability, even 
at energies where many partial waves contribute to the 
cross sections. The parameters can be determined exper- 
imentally, without having to rely on accurate PESs. The 
theory further relates ultracold ion-atom interactions to 
interactions at much higher temperatures. It is the be- 
ginning of a much broader program aimed at connect- 
ing ultracold interactions and reactions of all types, to 
their behaviors at much higher energies and/or tempera- 
tures that are more relevant to astrophysics and everyday 
chemistry. 

For simplicity and for purposes of making connections 
with existing formulations and setting benchmarks for 
further theoretical developments, we adopt here the so- 
called elastic approximation (see, e.g., Refs. [13, [Hj]) 
which ignores the hyperfine and isotope effects. Go- 
ing beyond this approximation will require a multichan- 
nel quantum-defect theory (MQDT) formulation of near- 
resonant charge exchange, along the lines of Refs. [13; HH 
and similar in spirit to recent works of Idziaszek et al. 
0, Q. Such a theory will be better understood after 
the conceptual and numerical developments of this work, 
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which relies only on the single-channel QDT [8[]. The 
multichannel formulation will be presented in a future 
publication. 

The paper is organized as follows. In Sec. Ill A| we out- 
line the application of QDT to the resonant charge ex- 
change of the type of 1 S+ 2 S. Section Hi Bl discusses two 
equivalent three-parameter descriptions derived from the 
QDT formulation. In Sec. IIII| the three-parameter QDT 
description is thoroughly tested for the case of resonant 
charge exchange of 23 Na, both through comparison with 
previous results of Cote and Dalgarno [10], and by com- 
paring our own numerical results with the corresponding 
QDT parametrization results. We conclude in Sec. [TV] 
with further discussions on the implications of our results 
in the more general context of ion-atom interactions. 



II. QDT OF RESONANT CHARGE EXCHANGE 

A. General considerations for 1 S+ 2 S type of 
systems 

Consider the system of an atom and its ion, one in a X S 
state, one in a 2 S state. This type covers resonant charge 
exchange of both Group IA (alkali), Group IIA (alkaline 
earth), and helium atoms in their ground states. Such 
a system correlates to two Born-Oppcnheimer molecu- 
lar states, characterized by 2 S+ U . The other, energeti- 
cally higher, electronic states can be ignored for collision 
energies much smaller than the first electronic excita- 
tion energy [TJj]. In the case of exact resonance or in 
the elastic approximation which ignores hyperfine struc- 
tures and/or the isotope shifts, the two coupled radial 
Schrodinger equations in the atomic basis become decou- 
pled in the molecular basis, reducing the understanding 
of resonant charge exchange to two single-channel equa- 
tions for t he g erade, q, and the ungerade, u, states, re- 
spectively da H, El 



h 2 d 2 h 2 i(i + i) 

~%i ZR 2 + 2[iR 2 



■ 9 ' U (R) = 



^x( e ) = ^E( 2/ + 1 ) sin2 (' 5 f-^)- (2) 

1=0 

It contains the physical concept that resonant charge ex- 
change is due to the phase difference between the gerade 
and the ungerade molecular states. For elastic and total 



cross sections, it is convenient to first define two single- 
channel "molecular" cross sections 



oo 



1=0 



l)sm 2 (Sn 



(3) 



in terms of which the total cross section is given by cr to t = 
(a 9 + <r u )/2, and the elastic cross section is given by 

0"e] = Ctot - 0"cx 

For a 1 S+ 2 S type of system with either a 1 S or a 2 S 
atom in its ground state, the potentials V 9jU (r) in Eq. ([1]) 
have the same leading term — C4/-R 4 at long range, where 
C4 > is given in atomic units by C4 = ai/2 with a± 
being the static dipole polarizability of the atom. Appli- 
cation of the single channel QDT for — 1/ R 4 type of po- 
tential Q gives an efficient characterization of the phase 
shifts Sf' u , leading to an efficient characterization of the 
resonant charge exchange process. 

Specifically, for a potential with a long-range behavior 
of V ~ —C n /R n , the tangent of the phase shift is given 
in QDT by 30] 



tan Si 



7 c(n) 



,c(n) 
fs 



Z c(n) K c 

gs 



Here K c (e, I) is a dimensionless short-range K c matrix 
[30j . Zx^\e s ,l) are universal QDT functions for —1/R n 
type of potential. They are functions of the angular 
momentum I and a scaled energy e s ~ e/sE, where 
s E = {h 2 /2^){1/ (3 n ) 2 is the energy scale and f3 n = 
{2^C n /h 2 ) 1 /<- n - 2 ^ is the length scale for the -C n /R n po- 
tential. Explicit expressions for Z^y~ 4 \ applicable to the 
polarization potential, are given in the Appendix. As ex- 
plained in Ref. 30] , Eq. (Q} includes not only the effect of 
long-range phase shift, but also effects of quantum reflec- 
tion and tunneling, which are the key differences between 
long-range potentials with n > 2 and those with n < 2. 

All cross sections for resonant charge exchange can be 
written explicitly in terms of tan^ ; 9 and tan 5™ for the 
gerade and the ungerade states. 



(1) o-ex(e) 
Here \i is the reduced mass, e is the energy in the center- 
of-mass frame, V g . u (R) represent the two potential en- anc j 
ergy curves for the gerade and the ungerade states, and 
u 9 t l u {R) are the corresponding radial wave functions for 
the Zth partial wave. 

In terms of the two phase shifts, S 9 ' u , for the gerade 
and ungerade states in partial wave I, as determined from 
the solutions of Eq. (p}, the total charge exchange cross 
section <r ox can be written as [13, [H, [HI 



1=0 



1) 



(tan<5f-tan^) 2 
(1 +tan 2 <5f)(l + tan 2 5?) 



„,. = !£ (2i+1) _'» 2 « 



1=0 



tan 2 (5f' u ) 



(5) 



(6) 



For sufficiently large I and away from a shape resonance, 
tan Si is independent of the short-range parameter and is 
given, for both q an d u states, by the Born approximation 
(see, e.g., Ref. Q) 



tan<5; 



(2Z + 3)(2Z + 1)(2Z-1) 



(7) 



This I /I 3 type of behavior for large I ensures convergence 
in summations over I in all total cross section calcula- 
tions. 
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The application of QDT allows the description of res- 
onant charge exchange in terms of the C4 coefficient, 
equivalently the atomic polarizability a±, and two short- 
range K matrices, Kg(e,l) for the gerade state and 
K°(e, I) for the ungerade state. Such a description is 
exact if the energy and the partial wave dependences of 
the K c s are fully accounted for. More importantly, QDT 
allows for efficient parametrizations of resonant charge 
exchange by taking advantage of the fact that the short- 
range K c matrices depend not only weakly on energy, but 
also weakly on the partial wave I for both atom-atom and 
ion-atom interactions [2{| |3(| • Through an example of 
Na + +Na, we show that even the simplest parametriza- 
tion, corresponding to ignoring the e and I dependences 
of the K c s completely, provides an accurate description 
of resonant charge exchange over a wide range of ener- 
gies, including energies at which tens of partial waves 
contribute. 



same set of parameters can in fact describe ion-atom in- 
teraction over a much wider range of energies, of the or- 
der of 10 5 S£;, as tested in the next section for 23 Na and 
expected to be qualitatively the same for all alkali-metal 
atoms. 

The two equivalent parametrizations are complemen- 
tary in terms of the physical understanding that they pro- 
vide. The parametrization using Kg, K°, and ai gives a 
more direct insight as to why it works over a wide range 
of energies. It is because Kg and K^ are short-range pa- 
rameters that are both insensitive to e and / [8, 29, 30]. 
The parametrization using a g i—o, a„; = o, and a± enforces 
the concept that the understanding of ultracold interac- 
tion immediately provides understanding of interactions 
over a much wider range of energies through QDT. This 
is because that embedded in the knowledge of the scat- 
tering lengths, a g i=o and a u i = o, are the knowledge of the 
Kg and K£ parameters, through Eq. {5}. 



B. Three-parameter QDT descriptions 

A three-parameter parametrization of resonant charge 
exchange for a 1 S+ 2 S system results from ignoring both 
the energy and the partial wave dependences of Kg(e,l) 
and K°(e, I), Specifically, it corresponds to the approx- 
imation of K°(e,l) « K c g {e = 0, 1 = 0) and K°(e,l) w 
K°(e = 0,1 = 0). Using Kg and K£ as the shorthand 
notation for the resulting constant K c s, we have one of 
the three-parameter parametrizations for resonant charge 
exchange, with two short-range parameters Kg and K°, 
characterizing the short-range ion-atom interaction, and 
one long-range parameter, C4 or the atomic polarizability 
cti, characterizing the strength of the long-range interac- 
tion. 

A mathematically equivalent three-parameter 
parametrization is in terms of two s wave scatter- 
ing lengths, a g i = o and a u i=o, for the gerade and the 
ungerade state, respectively, and the atomic polarizabil- 
ity a\. This is derived from the first parametrization 
by noting that the K c (e — 0,1 — 0), for both g and 
u states, are related rigorously to the corresponding s 
wave scattering lengths by [lH, [34| 



IR - f^b T 0--b) \ K c (0,0)+tan(^/2) 
a l=0 /p n - r(i + b) j Rc ^ Q) _ tan(7r6/2) 

where b = l/(n — 2). It reduces to, for n = 4, 

K c {0,0) + 1 



a;=o//?4 



K c (0,0) - 1 ' 



(8) 



(9) 



We note that in the context of the effective range the- 
ory [3514381 ] . the three parameters, a g i—o, a u i = o, and ot\, 
can only be expected to describe ion-atom interaction 
in the ultracold regime as characterized by e <C se, in 
which only the s wave makes a significant contribution. 
The QDT for ion-atom interaction asserts that the very 



III. THE EXAMPLE OF Na + + Na 

Low energy Na + +Na charge exchange for 23 Na has 
been studied in detail by Cote and Dalgarno in Ref. [l(| > 
within the elastic approximation. It serves as a proto- 
typical system to test the QDT formulation for resonant 
charge exchange, in particular the range of validity of the 
three-parameter description. 

In Sec. IIII A[ we make a preliminary evaluation of 
the QDT description by showing, visually, that a three- 
parameter QDT parametrization, using parameters as 
given in Ref. [10(, reproduces the cross sections of 
Ref. [Tpj ] . including all the resonance structures, with- 
out any knowledge of the short-range potential. A more 
detailed comparison is not possible since Ref. [10( made 
use of unpublished potential energy results by Magnier 
et al. that are unavailable to us. 

For a more detailed comparison between fully quantum 
numerical calculations and three-parameter QDT results, 
we construct in Sec. IIII Bl the 2 S+ and 2 S+ potential 
curves for 23 NaJ using the same procedure as prescribed 
in Ref. except by using the later published results 
of Magnier et al. [3§| in an intermediate region. These 
potentials are meant to be as close to those of Ref. as 
possible. Fully quantum numerical calculations of cross 
sections are carried out with these potentials and com- 
pared to the results of corresponding three-parameter 
QDT descriptions, and also to the earlier results of Cote 
and Dalgarno ^o| . 

In both sets of calculations, we take the sodium static 
dipole polarizability to be a.\ = 162.7 a. u. 5o| . the same 
as that adopted by Cote and Dalgarno [10| . The corre- 
sponding length scale for 23 NaJ is £4 = {2^Ci/h 2 ) 1 / 2 = 
1846 ao, where do is the Bohr radius, and the correspond- 
ing energy scale is seI^b = 2.21 fiK or Sfi/h = 46.05 
kHz. 
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Energy (a.u.) 



FIG. 1. (Color online) The total and the partial "molecu- 
lar" cross sections for the 2 Ej" state of 23 Na2~ from a three- 
parameter QDT calculation using a 3 ;=n = 763. 3ao, aui=o = 
7721.4a , and cti = 162.7 a.u.,all from Ref. [if}. It is to be 
compared with Fig. 2 of Ref. [lOj ] . The total cross section is 
obtained by summing over all partial waves until convergence. 
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FIG. 2. (Color online) The same as Fig. [TJ except that it is for 
the 2 Sj state of 23 Na+. This fi gure is to be compared with 
Fig. 3 of Ref. 



These s wave scattering leng ths, plus the Na polarizabil- 
ity of a\ = 162.7 a.u. [40| . give us all the parameters 
required for a three-parameter QDT description of reso- 
nant charge exchange, from which all relevant cross sec- 
tions can be calculated, without detailed knowledge of 
the potentials. 

Specifically, from the s wave scattering lengths of 
Eq. ([TTJ]) , we first calculate, using Eq. © , the short range 
K c parameters Kg and K^, and obtain 

Kg = -2.4095 , (11a) 
K c u = 1.6286 . (lib) 

In the three-parameter QDT description, they are taken 
as constants applicable at all energies and for all par- 
tial waves. These parameters, together with the QDT 
equation for the phase shift, Eq. (gj), give us all phase 
shifts and all cross sections. The results for the total and 
partial "molecular" cross sections for the 2 £+ and the 
2 E+ states are illustrated in Figs. Q] and [2] respectively. 
They are visually nearly identical to the corresponding 
results shown in Figs. 2 and 3 of Ref. Q3. All detailed 
features of the cross sections that are visible on the fig- 
ures are found to be at the right places judged by vi- 
sual examination. The results show, at least tentatively, 
that the three-parameter QDT description can provide 
an accurate account of resonant charge exchange over a 
wide range of energies, including all the complex struc- 
tures which in this case are shape resonances from a wide 
range of partial waves. 



B. More detailed comparison of QDT and 
numerical results 

For a more detailed comparison between fully quantum 
numerical results and three-parameter QDT calculations, 
we construct here a version of the 2 £+ and 2 S+ potential 
curves for NaJ. Numerical results are calculated using 
these potentials and compared to the QDT results corre- 
sponding to the same potentials. 



A. Comparison of QDT results with previous 
results 

Numerical results of the "molecular" cross sections for 
2 E+ and 2 £+ states, as defined by Eq. (|3]), are given for 
energies ranging from 10~ 16 a.u. to 1 a.u. in Ref. (Tol ) . 
The reference also gives the zero energy s wave scattering 
lengths for 2 E+ and 2 E+ states, 

a gl=Q = 763. 3a , (10a) 
a u i =0 = 7721.4oo . (10b) 



1. The potential energy curves adopted 

For both 2 T,f and 2 Slt" states of Nat , we use the ab 

i — i 

initio data of Magnier et al. [39J ranging from 5.0 ao 
to 20.0 a . Outside of this region, the potentials are 
extended using the same procedure as prescribed in 
Ref. [To| . in the hope of getting potentials as close to 
those of Ref. [l(| as possible. Specifically, for distance 
larger than 22.0 ao, we extended the potential by the 
asymptotic form of [10| 

Vg, u (R) = V disp (R) T V cxch (R) , (12) 
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FIG. 3. (Color online) Potential energy curves for 2 E+ and 
2 Ej states of Na^~ adopted in our numerical calculations. It 
can be compared to Fig. 1 of Ref. [lol ]. 



with =F for 2 S+ and 2 S+, respectivel y. The dispersion 



term and exchange term are given by (lCf 



V disp (R) 



i? 4 



V exch (R) = 2 ARae 



i? 6 

-bit 



B 

1+ R 



(13) 
(14) 



Except for different notations [4l|, all coefficients are 
taken to be same as in Ref. [Toj . which, in atomic units, 
are given by C4 = ai/2 = 81.35, Ce = 936.5, C$ — 
27069.5, A = 0.111, a = 2.254, b = 0.615, and B = 0.494. 
Potential energies between 20.0oo and 22.0a are inter- 
polated using a cubic spline [42| to make a smooth con- 
nection between the ab initio data and the long-range 
behavior. At short distances (R < 5.0a ), we extended 
the potential with an exponential wall as in Ref. [10] 



V(R) = We^q>{-wR) 



(15) 



with 



W = V(R)exp(wR)\ 



dlnV(R) 



5.0ao 5 



dR 



5.0a 



The resulting 2 E+ and 2 E+ potentials for Na^~, thus con- 
structed, are illustrated in Fig. [31 They have the same 
long-range behavior as those of Ref. [Ty] , and differ only 
slightly in the short range due to slightly different ab 
initio data adopted. 



2. Comparison of QDT and numerical results 

For the QDT calculations, we first calculate the param- 
eters Kg and K^, more specifically the Kg(e = 0, 1 = 0) 
and K^(e = 0, 1 = 0) from the potentials. The radial 
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FIG. 4. (Color online) Comparison of the charge exchange 
cross sections of 23 Na obtained from a three-parameter QDT 
description (dashed line) and from numerical calculations 
(solid line). 
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FIG. 5. (Color online) Comparison of the total cross sections 
of 23 Na + + 23 Na obtained from a three-parameter QDT de- 
scription (dashed line) and from numerical calculations (solid 
line) . 



wave function is matched to 

« d (r) = Mf: gl (r s ) - K c (e,l)g c eJ (r s )} , (17) 

at progressively larger R until the resulting K c converges 
to a constant to a desired accuracy [2^, [30, ■ Here f c 
and g c are the zero-ener gy Q DT reference functions for 
the — 1/R 4 potential [30l. l34j. We obtain 



K c 

9 



-1.5953 
0.25416 . 



(18a) 
(18b) 
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FIG. 6. (Color online) Comparison of charge exchange cross 
sections of 23 Na for two potentials both using three-parameter 
QDT description. Solid line: results for the potential of 
Ref. [ifjj] . Dashed line: results for our potential. 



From the K c s, the s wave scattering lengths can be ob- 
tained by substitution into Eq. ([9]). We obtain 

a g i =a = 423.51a , (19a) 
a M ;=o = -3104.8a . (19b) 

We note in passing that this method of calculating the 
scattering length converges at much smaller R and pro- 
vides more accurate result than by matching to the free- 
particle solutions or by matching the phase shift to the ef- 
fective range expansion [Io|, IH[ , especially for cases with 
ai =0 > Pa. 

Our numerical calculations of the phase shifts and cross 
sections are carried out using a log-derivative method 
[43l [iij . Figure [4] shows the comparison of the charge 
exchange cross sections obtained from numerical calcu- 
lations and from three-parameter QDT description. Fig- 
ure [5] shows a similar comparison of the total cross sec- 
tions. They both show that the QDT description and the 
numerical results are in excellent agreement for energies 
below 0.2 mK. For energy between 0.2 mK and 0.1 K, the 
QDT prediction still works well with the only discernible 
differences being due to shape resonances in high partial 
waves. A better QDT description of such resonances is 
possible, but is beyond the scope of this work that focuses 
on the simplest parametrization. Overall, the QDT pre- 
diction is satisfactory below 0.1 K, or about 10 5 S£. For 
energy higher than 0.1 K, the discrepancy between the 
two results grows larger. 

The considerable differences between the scattering 
lengths for our potentials, as given by Eq. (fl"9|) . and for 
the potentials of Ref. [l(| , as given in Eq. (fTUj) , are illus- 
trations of the sensitive dependence of cold or ultracold 
ion-atom interactions on the short-range potential. The 



FIG. 7. (Color online) Comparison of total cross sections of 
23 Na + + 23 Na obtained from three-parameter QDT descrip- 
tions using parameters corresponding to the potential of 
Ref. [l(| (solid line) and using parameters corresponding to 
our potential (dashed line). 



two sets of potentials have the same long-range behav- 
iors as characterized by Eqs. p^|) - (fl4l . and differ only 
slightly in the short range, to a degree that is visually 
indistinguishable on the scale of Fig. [3] Figures [6] and [7j 
which compare QDT results for the two sets of potentials, 
give a more complete picture of this dependence. They 
show that the interactions in the ultracold regime are 
the most sensitive to the short-range potential. Beyond 
the ultracold regime of e ^< s^;, the shape resonance 
positions remain sensitive to the potential and the QDT 
parameters over a considerable wider range of energies, 
of the order of 1 mK or about 1000s£. While the sensi- 
tive dependence of ion-atom interaction on the potential 
gradually diminishes at higher energies for the total cross 
section, it remains to a considerable degree for the charge 
exchange cross section. 

Fortunately, the sensitive dependence of ion-atom in- 
teraction on the PESs is fully encapsulated in a few (two) 
QDT parameters, as shown in Figs. [4] and [5] Instead of 
from the PESs, this small number of parameters can be 
determined from a few experimental data points, such as 
two resonance positions for Na + +Na, or two binding en- 
ergy measurements of vibrationally highly excited NaJ. 
Such a determination, in a similar manner as illustrated 
earlier for atom-atom interaction [29j, |45| , is further fa- 
cilitated by the concept of universal spectrum including 
the concept of universal resonance spectrum introduced 
in Ref. H for the — 1/R 4 potential. We look forward to 
the availability of experimental data to demonstrate such 
applications for ion-atom systems. 
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IV. DISCUSSIONS AND CONCLUSIONS 

In conclusion, we have shown that resonant charge ex- 
change of the type of 1 S+ 2 S can be accurately described 
over a wide range of energies using only three parame- 
ters, which can either be two short-range K° matrices, 
Kg and K°, and the atomic polarizability a\, or two s 
wave scattering lengths, a g i = o and a u ;=o, and ax- Since 
the polarizability is well known for most atoms, this is ef- 
fectively a two-parameter description. Everything else is 
described by analytic QDT functions for the — 1/i? 4 po- 
larization potential (see, Ref. @ and the Appendix). In 
the case of 23 Na, excellent agreement between the QDT 
parametrization and numerical results is found from K 
all the way through 0.1 K, including all resonances within 
this range. To put this energy range into perspective, we 
note again that e/kB — 0.1 K corresponds roughly to 
e s = c/se ~ 10 5 . At this energy, one can estimate that 

there are at least \f2tj 4 ~ 25 partial waves contributing 
to the cross sections. Such a simple parametrization over 
such a wide range of energies is made possible in QDT 
not only by the energy insensitivity of the short-range pa- 
rameters, Kj: and K^ ll but also by their partial- wave in- 
sensitivity [g, H^, H3| • The energy insensitivity is ensured 
here by the large length scale separation as is typical for 
ion-atom interactions. More quantitatively, it is reflected 
in /3 6 /^4 « 6.56 x 10~ 3 for Na, where (3 6 = (2fiC 6 /h 2 ) 1/4 
is the length scale associated with the —Cq/R 6 term of 
the potential in Eq. (p~3|) . This value gives an order-of- 
magnitude measure of the length scale separation that is 
representative of all alkali-metal atoms. The partial- wave 
insensitivity is ensured by the combination of length scale 
separation and the smallness of the electron-to-nucleus 
mass ratio [1^, [34[ . 

The example of 23 Na has also served to illustrate the 
sensitive dependence of cold or ultracold ion-atom inter- 
actions on the PESs. While one can always construct 
the potentials for ion-atom systems, their accuracies are 
generally far from sufficient in predicting cold collisions, 
and the related vibrationally highly excited molecular- 
ion spectrum [|. The QDT deals with this difficulty by 
encapsulating this sensitive dependence into a few short- 
range parameters that can be determined experimentally. 
The simplicity of the resulting description has the follow- 
ing implications, (a) Since there are only a few parame- 
ters, they can be determined from very few experimental 
data points, (b) Since the parametrization works over 
a wide range of energies, it allows the determination of 
the parameters from measurements of structures (either 
resonance or binding energy) away from the threshold, 
where they are much further separated [8[ and can be 
resolved with much less stringent requirement on cither 
the energy resolution or the temperature [20]. This can 
be an important consideration for ion-atom interactions, 
where going below millikelvin has been proven to be dif- 
ficult experimentally, (c) Such a parametrization offers a 
systematic understanding of a class of systems. For ex- 
ample, the parametrization for Na + +Na works the same 



for all resonant charge exchange processes of the type of 
1 S+ 2 S. Different systems differ only in scaling as deter- 
mined by the atomic polarizability a.\, and the two short- 
range parameters, (d) The parameters that are used to 
characterize interaction in the absence of any external 
field also characterize the interaction in the presence of 
external fields, thus rela ting field-free interactions to in- 
teractions within a field [3l|, |46J . 

This work represents only the simplest QDT descrip- 
tion for resonant charge exchange, with a goal of estab- 
lishing key qualitative features that form the conceptual 
foundation for further theoretical development. Future 
works will include rigorous treatments of the hypcrfine 
and isotopes effects, more accurate treatment of reso- 
nances in high partial waves which give rise to only de- 
viations of any significance in Figs. U and El and further 
extension of the range of energies to room temperatures 
and above. A realization of these goals, with a theory of a 
few parameters, will greatly facilitate the incorporation 
of accurate ion-atom interaction data into simulations 
of not only cold plasmas @, but also plasma systems 
of interest in astrophysics and everyday chemistry and 
technology. 
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Appendix: The Z c{n) matrix for potential 

The Z c ( n ) matrix is well defined for any —1/R n type of 
potentials with n > 2 [3(| • For the polarization potential 
corresponding to n = 4, its elements are given by 

7 c(4), n _ COSK^ - ^o)/ 2 ] 

*fs ^ l >~ 2M e .,cos(7ri//2) 

x {1 - (-l) l Mlt tan[7r(^ - i*)/2]} , (A.l) 
7 c(4), n _ cos[tt(;/ - z/q)/2] 
* f c - 2Af £s; cos(W2) 

x {tan[^ - vo)/2] - (-l)'Af 2 ,} , (A.2) 

7 c(4)/ n _ COs[7T(i/- i/p)/2] 

^ gs - 2M eJS m(nu/2) 

x {1 + (-l)'M 2 , tan[7r(t/ - vo)/2]} , (A.3) 
7 c{A) ( n _ cosK;/ - vo)/2\ 

^gc K^,i) ~ 2Mejsin ( nl//2 ) 

x {tan[^ - vo)/2] + (-l)'M 2 ,} . (A.4) 

Here Vq = 1 + 1/2, v is the characteristic exponent for the 
— 1/i? 4 potential (corresponding to the modified Mathieu 
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equation [Tm - 1491 ] ) . and M e j is one of its QDT functions. 
Their evaluations have been discussed in Ref. Q. To- 
gether with Eqs. (jA.ljl - (|A.4| ) for the Z c ^ matrix, they 
are all the requirements for the implementation of QDT 



for resonant charge exchange. The derivation of Z cl ^\ 
and other aspects of QDT for — 1/i? 4 potential, will be 
presented in a separate publication. 
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